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The theory developed in an earlier paper, for an accurate reference orbit of an artificial satellite, 
is first slightly modified, so as to prepare the way for a treatment of zonal harmonic perturbations. 
Delaunay variables are next introduced, by means of certain linear combinations of the action variables, 
along with their canonical conjugates. Application of the von Zeipel method then permits the cal- 
culation of the most important zonal harmonic perturbations. These arise from the third, with co- 
efficient y 3 , and the residual fourth, with coefficient J*+J\. The accuracy of the secular and short- 
periodic effects is through terms of order }\ and that of the long-periodic effects is through terms of 
order ]i. Since the reference orbit itself, with its exact secular terms, takes care of all but 0.5 per- 
cent of the deviation of the earth's gravitational field from spherical symmetry, the overall secular 
accuracy of the final orbit surpasses that of other second order theories. The results are compared 
with those of Kozai. 



1. Reference Orbit 
The author [Vinti 1959a, b] has introduced a potential 

r=-/xp(p 2 + c 2 7?V (1-00) 

that can represent accurately the gravitational field .of an oblate planet. Here p, is the product of 
the gravitational constant G and the planet's mass, c is an adjustable length, and p and r) are 
oblate spheroidal coordinates, defined by the equations 

X + iY= r cos 6 exp i<f> = [(p 2 + c 2 ) (1 - rf)]* exp i<f> (1.01) 

Z = r sin 6 = pr). (1.02) 

If an artificial satellite is at the field point, r, 0, and </> are respectively its planetocentric distance, 
declination, and right ascension, and A, F, and Z are its rectangular coordinates, OZ being along the 
planet's axis and OX pointing toward its vernal equinox. 
If r> is the equatorial radius, the true potential is 

V=-fjLr~ l [l -^(re/ry'JnPnis'm 0)] + tesseral harmonics, (1.03) 

n = 2 

where only the products r£/„ need to be known. That is, differences in the definition of r e , when 
noncircularity of the equator is taken into account, can be reconciled by small adjustments of the 
7's. Then if 

c = r e J*, (1.04) 
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V represents exactly the zeroth harmonic — /x/r and the second harmonic and also gives higher 
even harmonics, characterized by 

J r 2 m = (-l) m+ Vr. (1.05) 

In particular it gives ,/4+./ 2 — 0, as compared with observed values for the earth ranging from 
-(0.9)10- 6 to (0.4)10- 6 (Kaula 1962; King-Hele, Cook, and Rees 1963). Consequently it accounts 
for about 99.5 percent of the deviation of V from the value — /x/r corresponding to spherical sym- 
metry. It thus accounts almost completely for the flattening of the earth, leading to a geoid that 
never departs by more than about 30 m from the true sea-level surface. 

For the drag-free motion of an artificial satellite the potential (1.00) leads to a separable prob- 
lem, which has been worked out analytically [Vinti 1961 a,b, 1962]. This solution, holding for 
all angles of inclination and containing no critical inclination or long-periodic terms, gives secular 
terms exactly by means of rapidly converging infinite series and short-periodic terms correctly 
through order J\. We call this orbit corresponding to (1.00) the reference orbit. For such a ref- 
erence orbit error can never accumulate, because of the exactness of the secular terms, and the 
periodic terms can be in error only by amounts of the order 73, i.e., by about 1 part in 10 9 , since 
/ 2 = (1.08)10- 3 forthe earth. 

2. Zonal Harmonic Perturbations 

For a satellite of the earth, if its orbit is high enough so that drag is small and low enough so 
that the moon's effect is small, the above reference orbit ought to hold rather well for a good many 
revolutions. (I purposely choose vague words here, since numerical comparisons are still incom- 
plete.) Eventually, however, the actual orbit will deviate more and more from such a reference 
orbit, because of the neglected forces. These include forces arising from drag, meteoritic impact, 
radiation, electromagnetic fields, the sun, and the moon, and the neglected part of the earth's 
gravitational potential, corresponding to (1.03) minus (1.00). Since the expansion of (1.00) in zonal 
harmonics is 

V =-tMr-^(r e /rr m (-J 2 rP 2m (sm 0), (2.00) 

m= 1 

this difference is 



V-V' = ixr- 1 



jfjiP&in d) + (^\j 4 + jl)P4(sm 6) + fej/^rfsiii 0) 



+ ^W 6 -4)P6(sin0)+. . 



+ tesseral harmonics. (2.01) 



Of these forces the most important, for any satellite with a large ratio of mass to area, are the forces 
corresponding to J 3 and J 4 + J\ in (2.01) and drag, which as determined empirically may include 
effects of meteoritic impact. For a double satellite [Langer and Vinti 1963] only (2.01) and the 
lunar-solar perturbation remain. 

The purpose of the present paper is to devise a method for correcting for the effects of any of 
the zonal harmonics in (2.01). The first example considered is the residual fourth harmonic, with 
coefficient J4 + J\. This harmonic leads not only to short-periodic effects and secular effects, but 
also to long-periodic effects depending on a resonance denominator 1-5 cos 2 /, giving rise to a 
critical inclination 7 = 63.4°. The second example considered is the third harmonic, with coeffi- 
cient 7 3 . This gives rise only to short-periodic effects and to long-periodic effects without singu- 
larities, so that it is qualitatively less interesting. Because of its greater magnitude, however, 
J 3 being about (— 2.4)10" 6 and I/4 + /H being probably somewhat less than (0.5)10~ 6 [Kaula 1962; 
King-Hele, Cook, and Rees 1963], it leads to somewhat larger periodic effects. 
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3. The Dynamical Problem 

Our problem is thus to find the motion of a satellite, taken to be of unit mass, when the Ham- 
iltonian is 

F = -7>^p(p 2 + cV) _1 +^i, (3.01) 

where T is its kinetic energy and where, for the residual fourth harmonic 

Fi = -/xr|r- 5 o- 4 P 4 (sin 6) (3.02) 

(H=Ja + JI (3.03) 

(We have here reversed the sign of the Hamiltonian, to agree with the usual practice with Delau- 
nay variables.) In carrying out this solution we shall use the results and notation of the solution 
[Vinti 1961 a,b, 1962] for the reference orbit, for which F' =0. If in (3.02) we then put 

r=a(l-e cos E) = a{\ -e 2 )(l + e cos v)~ l (3.04) 

sin 6 = sin / sin(i; + j8 2 ), (3.05) 

the expression (3.02) for Fi will be correct through order J\. To this order of accuracy the anom- 
alies E and v may be given by the quasi-elliptic expressions 

E~M s + Eo (3.06) v~M s + vo. (3.07) 

Note that (3.04) corresponds to r~ p and (3.05) to other approximations of zeroth order in y 2 , viz, 

ili~ili 8 + *lto~Ms + p2 + vo~v + P2 (3.08) rj-sintf. (3.09) 

Such an order of accuracy will result in errors of order J\ for those secular and short-periodic 
effects which are produced by the perturbing potential (3.02) and of order J\ for the correspond- 
ing Long-periodic effects. This perturbation (3.02) represents about 0.1 percent of the departure of 
the earth from sphericity. The solution for the other harmonics in (2.01) will have the same 
accuracy. However, since all of these higher harmonics represent only about 0.5 percent of the 
earth's departure from sphericity, their lower accuracy, as compared with that of the reference 
orbit which has already accounted for 99.5 percent of this departure, should not result in serious 
cumulative errors. 

In doing the perturbation theory, the first canonical variables that come to mind are the 
Jacobi "constants," viz, the a's and jS's of the reference orbit. When the reference orbit is ellip- 
tic, however, their shortcomings are well known and they lead to the same troubles in the present 
problem, giving rise to Poisson terms, linear in the time, in the variations in a.\ and a 2 . 

The next set of canonical variables that one might try is the set generated from the a's and /3's 
by the generating function 

S'=-ai* + M-2a0-*j8J + a^ + atf;. (3.10) 

If we define no by 

fjL=nlal a = — J/x/«i, (3.11) 
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the resulting canonical variables are 

L = (fjLa )* / = rco(*+/3i) 

a 2 j8 2 (3.12) 

canonical with respect to the Hamiltonian 

F=ML*+Fi. (3.13) 

When the reference orbit is elliptic, this set is the same as the fast Delaunay set [Garfinkel I960]. 
One may then attempt to apply the von Zeipel method in the way successfully used by Brouwer 
[1959] and Garfinkel [1959], first eliminating short-periodic terms and then proceeding to eliminate 
long-periodic terms. One finds, however, that the corresponding generating function S*, which 
ought to be of the first order in the parameter cr 4 = J\+ J\, must then satisfy 

dSf/djS; = zeroth order in a 4 . (3.14) 

One may alternatively eliminate short-periodic and long-periodic terms simultaneously, but 
one then obtains a Poisson term of the form v' sin 2/3' 2 in a 2 — a! 2 . Since v' has a secular part, such 
a result would appear absurd, since the "constant" a 2 , which ought to have only a small periodic 
variation, would then increase indefinitely with time. 

These difficulties are examples of the failure of the von Zeipel method whenever the following 
conditions both hold: (1) the perturbing potential has a long-periodic part of the first order in the 
perturbation parameter, and (2) the canonical variables are such that the unperturbed Hamiltonian 
depends only on L. 

To obtain a successful set of variables, we may proceed as follows. Let qf and p?, i= 1,2,3, 
be the coordinates and momenta p, 77, </>, p p , p^, p<t> corresponding to the unperturbed problem (the 
reference orbit), with Hamiltonian F = Fq. Also let 77, wf, i= 1,2,3, be the corresponding action 
and angle variables. 
Then 

P?=aS(&A&>l,&#/dd (3.15) 

«#=dS(& qt ql&jlMVdJl (i=l,2,3) (3.16) 

where S is the Hamilton-Jacobi function of the unperturbed problem [Vinti 1959b], with the Jacobi 
a's replaced by the fl. Here 

fl = jrtdtf=jl(aua* 9 ati. (3.17) 

Now let qt, pi, i = 1,2,3, be the coordinates and momenta corresponding to the perturbed problem, 
with Hamiltonian F = Fo + Fi. Introduce new variables 7'/, wu i = 1,2,3, by means of the canonical 
transformation 

Pi = dS(q u <?2, q.u ji,J2,jz)ldqi (3.18) 

Wi = dS(q u q 2 , qzjujtj*)ldj u (i= 1,2,3) (3.19) 

where S is the same function of the q\ and j) that the above Hamilton-Jacobi function is of the 
g? and;^. Then the w\ and j\ are canonical variables, satisfying the equations 

djildt = dFldwi^O (3.20) 

dwildt = - dF/dji 7^ constant (3.21) 
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They are thus not action and angle variables, since the ji are not constant and the w% are not linear 
functions of t. Moreover 



ji * Jpidqi, (3.22) 

in contradistinction with (3.17). 

It pays to go further, however, and introduce still another set of variables, a new Delaunay set 
L, G, //, /, g, h, by the transformations 

27rL =71 +72 +73 s S n <*3 1 = 2ttw\ 

2ttG=j 2 +73 sgn a 3 (3.23) g= 2tt(w 2 - w x ) (3.24) 

2irH= z js h = 2tt{w>>, — w> sgn #3), 

where sgn Q£3=±l respectively for a direct orbit or a retrograde orbit. To verify that they are 
canonical, note that 

Ldl + Gdg-\- Hdh =j\dw\ -\-j2dw2 +jadwz. (3.25) 

They were used by Izsak [1962] in his application of the author's reference orbit to the problem 
of the critical inclination. From (3.23) we now have 

7, = 2tt{L - G) 72 = 2tt(G - H sgn a 3 ) 73 = 2ttH. (3.26) 

4. The New Delaunay Set 

The functional relations among any of the quantities a ? -, 7'/, w% for the perturbed problem are the 
same as those connecting cuf, 7*?, uf\ for the unperturbed problem. We may therefore usually drop 
superscript zeros and depend on the context for the meanings of the quantities. 

From the author's paper fVinti 1959 b] we now find 

7i = 2 \ f,2 ^dp = 2Tr\ti(-2a l )-* - a 2 ] + 0(J 2 ) (4.00) 

J pi d P 

72 = 4 T d V = 27r(a 2 - a 3 sgn a 3 ) + 0(J>) (4.01) 

J Q Of) 



Next put 



Then, since 



find 



73 = 2ttoj3. (4.02) 

a r8 = dctrldjs (4.03) v { = doPJdfi (4.04) j n s dj r /da s . (4.05) 

if? ^=1^-, (£=1,2,3) (4.06) 

A^ hi i (i(Yi- firth 



i = l 



dji da k dak 



^- = * + j3i =7hm;i +721^2 (4.07) 

d(X\ 

— = j8 2 =712^1 +722^2 (4.08) 

dc*2 

T~ = j8 : , =713^1 + /23W2 + 27Tii;:r (4.09) 

5 0:3 
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With the aid of (3.24) these equations become 

2ir(t + p 1 )=j ll l+j il (l + g) (4.10) 

27rp 2 =j 12 l+J22(l + g) (4.11) 

2tt/3, =j v ,l + (/ 23 + 2tt sgn a.,)(/ + g) + 2ttA. (4. 12) 

The quantities y,s occurring here are given explicitly as functions of the as in eqs (7.16) through 
(7.21) of an earlier paper [Vinti 1961a]. 

The constant orbital elements in the perturbed problem are then the constant parts a", e", 
t?o of a,e,7? , along with the initial values JJ, gj, AJ of the secular parts of /, g, h. The corresponding 
Hamiltonian F is given by 

F = F (L,G,H) + F l (4.13) F = -a l , (4.14) 

so that 

L = dF/dl l=-dF/dL 

G = dF/dg (4.15) g=-dF/dG (4.16) 

// = dF/d/i = 0. h = - 8F/dH. 

With the use of (3.23), (3.26), and (4.04), we find for the unperturbed problem that 

dF da\ JL da, dji 

L- dfo_ Ha, _ ^ da, dji 
h -~~dH-~dH-& W 3H = 2Mvs ~ Vl Sgn as) 

and thus for the unperturbed problem that 

/ = /()+ 2lTV\t 

g=go + 27r(v 2 -v l )t (4.18) 

h = ho + 27r(^3 — ^2 sgn a 3 )i, 

where we have dropped the double primes from / , go, ho. 

Before finding the effects of the perturbing potential, it is desirable to change the algorithm for 
the unperturbed problem, given in an earlier paper [Vinti 1961 a, pp. 197-200], so that the constant 
orbital elements become a, e, /, /o, go, and {$3. {@3 is better than /i , as we shall see later.) To do 
so, insert (4.10), (4.11), and (4.18) into eqs (8.2) and (8.3) of that paper and carry out the same 
process that was carried out there. One finds 

M 8 = lo + 2iTPit (4.19) 

ipM = k+gQ + 2irv2t. (4.20) 

With these new definitions of M x and \\) s the algorithm then becomes the same as in the earlier 
paper. Note, however, that the restriction on the angle of inclination / has been removed 
[Vinti 1962]. 

5. Variations in the Spheroidal Coordinates 

It is convenient to derive here the variations in p, 77, and <f> that will arise from the variations 
produced by the perturbing potential in the Delaunay variables. From the later solution of the 
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canonical equations (4.15) and (4.16) we shall find that the rath zonal harmonic produces variations 
in the Delaunay variables of the forms 

8/. = (TmLm 81 = (Tmlm + CTmJ^lm 

8G = (T m Gm + (TmJ^Gm (5.00) 8g = dmgm + (TnJ^gm (5.01) 

8H = 0. 8h = (T m hm + CTnJ^hm. 

Here 0"3=J3 and CT4 = J4 + J}- The terms o- w L m and (T m G m are short-periodic of order j\\ the 
products of (TmJl 1 with G m , / m , g m , and h m are long-periodic of order J 2 . For ra = 3 the products of 
C73 with / 3 , g3, and /i 3 are short-periodic of order J$. For ra = 4 

U = Ml "T" *42 

g 4 = g 41 +g 42 (5.02) 

^4 = ^41 + ^42, 

where the products of 0-4 with / 4 i, g4i, and A41 are short-periodic terms of order 7^ and where the 
products of o~4 with / 4 2, #42, and /142 are secular terms of order J\. 

From 8L and 8£ we can find the variations of the /s, then of the a's, and finally of the ele- 
ments a,e, and 170 — sin /. From these and from 8/, 8g, and 8h we can then find 8E, 8v, 81//, and 
8x and finally the coordinate variations dp, 8rj, and 8(f). To obtain the 8/'s, use (3.26) and (5.00). 
The results are 

871 = 2lT(cr m L m — (TrnQm ~ (TmOmlJl) 

8j z = 2rr{cr m G m + (T m G m IJ 2 ) (5.03) 

8/3 = 0. 

To find the 8a's use 

S«r=i^8/V (5.04) 

Within the accuracy of the calculation, the coefficients da r /djj are needed only through order 
Jl Thus by (5.03), (5.04), and (4.04) 

8a x = 2iTV\(T m L m + 2tt{v 2 — v\)cr m G m + 27r(^ 2 — V\)(T m G m IJ2- (5.05) 

Then, by (5.03), (5.04), and (4.03) 

Sa 2 = 2TTCL2\(TmL m + 27r(a:22 ~ OL2\)(T m Gm + 277(a 2 2 ~ OLi\)V m G m \ J <L* (5.06) 

Also, by (4.02) and (5.03) 

8o: 3 = 0. (5.07) 



To find a 2 ] and a>2 note that 



so that 



^ do2 djj_ __ 802 _ 
£i dji da 2 da 2 
3 da> djj = da 2 = Q 
£1 3fl d a i dc*i 



a2X/21 + a22722 = l 

a2i,/ii + 0:22712 = 0, 



(5.08) 



(5.09) 
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with the solution 

a 2 i = -\/2i/A Q!22=Jii/A, (5.10) 

where 

A =7*11722—712721. (5.11) 

From an earlier paper [Vinti 1961 a, p 189], it then follows that 

an = OC/2) « 22 = (277)" 1 + 0(A) (5.12) 

and also that 

27m = /i + Ot/ 2 ) 

(5.13) 

V2 — Vl=0{J 2 ), 

where 

n = ^a~^ 2 . (5.14) 

On carrying secular and short-periodic variations only through order J\ and long-periodic perturba- 
tions only through order J 2 , it then follows from (5.04) through (5.08) and (5.12) through (5.14) that 

dai — ncrmLm (5.15) 

da 2 = (T m Gm + <T m GmlJ2 (5.16) 

60:3 = 0. (5.17) 

In finding the variations of other quantities we may drop all terms in their defining equations of 
order /*, where k ^ 1. To show this, note that if P= 0(J%) and 8P = <j m P\ + <JmJ 2 x P\ where P\ and 
P\ are both of order /J, then 



F 2' 



8UiF)=J%<rmPi + o-rJ^Pi). (5.18) 

Thus if A: i= 1, /| o- m Pi is a secular plus short-periodic term of order A; + 2 i? 3 and ]\ (TmJ~ l Pi is a 
long-periodic term of order k-\- 1 i? 2. This proves the statement. 

To obtain the variations of a, e, and 170 we may thus use 

(5.19) 



JL 

2a x 


1- 


9 20:1^2 

— p^ 555 


2 

, 2 «3 
CK 2 



the neglected terms being of order / 2 . There follow 

da =— - 8ai = Sai 

2af /x 



e8e = 4s«i+^8a 2 (5.20) 



1706170 = —? 8a 2 - 



With use of (5.15), (5.16), (5.19), (5.20) and the relation 

a\ = fxp-\- OU2), we then find 



8a = 2VmLrn (5 21) 



.,(,)- 



DTI \ P ) ~ 

de = — cr m L m — (aey l \— (<T m G m + (TmJ^Gm) (5.22) 

fie fi 

Syo = VoKw)'* (1 — T7o)(°~mG m + CTmJ^Gm). (5.23) 
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The element a thus has only a short-periodic variation while the elements e and rjo = sin / have 
both short-periodic and long-periodic but no secular variations. 

To find 8E, 8v, and 81// we first insert (4.10) and (4.11) into eqs (8.2) and (8.3) of the earlier paper 
[Vinti 1961 a], rejecting all terms of order J%, where k ^ 1. We find 

E - e sin E = I + 0(/ 2 ) (5.24) 

x\, = v + g+0{J 2 ) (5.25) 

Equation (5.24) gives 

8E = (1 - e cos E)- X (8l + sin E8e). (5.26) 



To find 8v , we use the anomaly connection 



v /l + e\* E 



tan 2\F^/ tan 2' (5 ' 27) 

from which there follows 

= (1 - e 2 )* (1 - e cos £ , )- 1 [8£* + (1 - e 2 )~ l sin £6>]. (5.28) 

Then 

8ifr = 8v + 8g. (5.29) 

The variations of the spheroidal coordinates 

p = a(l-e cos E) (5.30) 

77 = 770 sin t// (5.31) 

are then 

8p = (l — e cos E)8a — a cosE8e + ae sin f E8E (5.32) 

677 = sin 1//8770 -1-770 cos 1//81// (5.33) 

To find the variation 8<fi of the right ascension, we note that by eq (8.49) of the earlier paper 
[Vinti 1961 a], 

<t> = fa + X sgn a 3 + OC/2), (5.34) 

where 

tan x = (l -*?!)* tani// (5.35) 

cos x = (1 "~ Vo sin 2 1//)~ % cos 1//. (5.36) 

Equation (5.36) was given in Vinti [1961 b]. Also by (4.12) and the results 

713 = OC/2) 713 4- 2tt sgn a 3 = 0(J 2 ) (5.37) 
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from page 189 of Vinti [1961 a], we find 

j3 3 = A + OC/ 2 ), (5.38) 

so that 

8(3 3 = 8h. (5.39) 

From (5.35) and (5.36) we find 

S X = (1 - Vl sin 2 i/O-'Kl - *?§)* ty- | (1 - t? 2 )-* Vo sin 2i//8t ?0 ]. (5.40) 

Then, from (5.34), (5.39), (5.40), and (5.23), it follows that 

80 = 8h + (1 - t? 2 sin 2 1//)" 1 cos /[Si// - § (/Lip)-* sin 2xjj(cr m G m + oW^C,,,)], (5.41) 

where we have used (1 — r)D* sgn a 3 = cos /. 

It is well to note here that /3 3 is just as useful an orbital element as h . 

6. Solution for the Perturbed Delaunay Variables 

We now have to solve the canonical equations (4.13) through (4.16), the perturbing potential 
being 



Ft =— /xre 4 r ^/Vsin 0) 



(6.00) 



in the case of the residual fourth harmonic. To obtain secular and short-periodic variations through 
order ]\ and long-periodic variations through order J>, it will suffice to use elliptic approximations in 
(6.00), since o- 4 = J a + J\ = 0(7 2 ). Thus in (6.00) we may put 



sin 6 = 7] = 7)o sin i// = 770 sin (v + g) 
r = p = a(l — ecosE) = a(l — e 2 )(l + ecos v)~ l 



(6.01) 
(6.02) 



and we may use 



Then 



Ft = 



- jAri(r 4 



8<z 5 



l = E — e sinE 
dl = (r/a) 2 (l-e 2 )-*dv 



3-15t)§ + 



105 



I- 1 +577 2 (3-^ 2 )P) cos2i// + 



35 Ja 

T^7. 



cos 4i// 



The part of F\ independent of / is then 

F im - (2.)- j>*-^ (W)-* /; [(3- 15„> + f «) g 

+ 5tjS(3-|t,s)(^ 3 cos 2^ + ^1,3^) 



(6.03) 
(6.04) 
(6.05) 



(6.06) 



cfo (6.07) 
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From (6.02) and (6.05) it follows that 



f(°)' 



^-277-(l-e 2 )- 3 ^H-|e 2 



r 

J It 



a\ :i 3 

-J cos 2\fj dv = -7re 2 (l — e 2 ) -3 cos 2g 



m 



cos 4i// dv = 0. 



Equations (6.07) through (6.10) then result in 
where the constant part is 



Fi c =- 



^(l-e 2 )-(l + |e 2 )(3-15, f) + f^ 

and the long-periodic part is 

Fu^-^g^d-er^vl^-lvl) cos 2g 

The short-periodic part is then given by 

Fu — F\— Fie — Fi p , 
so that 



/xr P o- 4 17., 1C o i 105 4 



5 / 3 

_ (1 _ e 2)-7/2/ 1+ | e 2 



From (3.23) and (4.00) through (4.02) we now obtain 

L = fi(-2ai)-* G = a 2 H=a :i , 
with neglect of terms of order J>>. To the same accuracy 

a = —%/JLa~ l al = /jia(l — e 2 ) a 2 a.7 2 = 1 — rj 2 . 
Thus, with neglect of terms of order J>, 

a = L 2 ijl- 1 l~e 2 = G 2 L~ 2 tj 2 = 1 - H 2 G\ 



(6.08) 
(6.09) 
(6.10) 

(6.11) 
(6.12) 

(6.13) 
(6.14) 



+ 5ij8(3~i|S) Q 5 cos2i^-|e 2 (l-e 2 )- 7 / 2 cos2^ + y t^)* cos 4i//j- (6.15) 



(6.16) 



(6.17) 



(6.18) 



as expected. When we later take derivatives of a generating function with respect to L, G, and H, 
we shall have to use the expressions (6.18) to replace them in the results by the elements a, e, 
and 770. 

In solving the canonical equations (4.15) and (4.16) we first make a canonical transformation 
to new canonical variables L', C, //', /', g' , and h', so that the new Hamiltonian F* shall be 
independent of /' and h' . This first step will yield the short-periodic effects. To carry it out, 
introduce the generating function 



S(L',G f ,H',l,g,h) = So + S l (L\G',H',l,g), 



(6.19) 
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where 

S* = L'l + G'g+H'h (6.20) 

and where Si is to be of the first order in cr 4 . Then follow the von Zeipel method, as applied by 
Brouwer (1959) to satellite orbits. 

On splitting F* into parts F$ and F*, we then have 

F (L, G, #)+Fi(L, G, #, /, g) = F$(L\ G', //') + Ff(L\ G\ H\g') (6.21) 

along with the following relations connecting Si and the primed and unprimed variables: 
j _ dS _ , dSi j, _ dS _ j . dSi 

^ dS ^. . dSi /ir ~~ x , dS . dSi _ __ 

G= ^ =G+ ^ ,6 - 22) *=ac^ + ^ < 6 - 23 > 

U—dS —II' Lf— §S _ L I ^Si 

Since the new Hamiltonian will not contain /' or h\ we shall have H' =H = constant and L' = con- 
stant. Insertion of (6.22) and (6.23) into (6.21) then gives a partial differential equation for Si: 

4' + f, G ' + f,H) + f ' 1 ( i ' + f, C ' + f.H./ + f„ + §) 

»FJ<L\C',fl)+Ff(z.',G',fl,«+J^Y (6.24) 

Taylor expansion in the neighborhood of L' , G', //, /, and g, with rejection of terms beyond the 
first order in cr 4 , then gives 

F <) (L',G',H) + ^^ + ^^+FiL',G',H,l,g) = F*(L',G',H} + Ft(L',G',H,g). (6.25) 

The zeroth order terms lead to 

F^L',G',ff) = F*(L',G',H) (6.26) 

and the first order terms to 

^%+^§± + FUL\G\m + F ip (L',G',H,g) + Fu(L',G',H,l,g) = Ff(L',G',H,g') 

(6.27) 

In writing down (6.27) we have used (6.14) to express Fi as a sum of the terms Fi c , Fi p , and Fu and 
we have replaced g by g' in Ff, a permissible step involving an error of order cr\. In (6.27) the 
terms independent of / then yield 

Ft = F 1C (L\ G', H) + F lp (L\ G\ //, g\ (6.28) 

so that (6.26) and (6.28) together provide the new Hamiltonian F*. The remaining terms, depending 
on /, then yield 

dFodSi dFodSi r , T , n , j, , 

ar^ + aG 7 ^ = " Fl((Z " G ' // ' / ^ ) ' (6 - 29) 
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a partial differential equation for S\. With use of (4.17) it becomes 

2ttv x -77 + 2tt(v 2 - v\ ) — 1 = Fi/(L', G', //, /, #) 



(6.30) 



Since v% — v\ is of order J 2 and since Si is to be of order cr 4 = J 4 + /|, it follows that 2tt(^2 — Vi)dS\ldg 
will be of order /f and is thus to be rejected. Thus 



27mS, = J FiKL', G\ //, /, g)d/ + 4>(g). 



(6.31) 



To evaluate the integral in (6.31), apply (6.15), (6.02), (6.04), and (6.18). Since Fu has a factor 
vi — OUl)* we can niake a number of approximations at this point and still achieve our desired 
accuracy. These are: drop the primes from L and G in (6.31), place 



27TVx = n = in*a- 3l \ 



(6.32) 



express a, e, and 770 by means of (6.18), drop the primes from L and G in calculating dSJdL' and 
dSJdG, and finally use (6.18) again to replace L, G, and H in the final formulas by a, e, and 7/0. 
We obtain 



- e 3 sin v + — e 2 sin 2i> + — sin 3t> 
4/4 12 



S t = (T4Q J (l + § e 2 ) (v - /) + (3e + 

+o- 4 (? 2 [|e 2 (2;-/)cos2^ + |-sin(t;--2^ + ^e+^ 3 



♦G-' 



-e 2 j sin (2t> 



( I + !) sin (3v + 2g) + H sin (<fc + 2tf) + ^ sin (5v + 2*)] 



+ 0-4^3 



3 Q / 3 

g- sin (i> + 4g) + - e 2 sin (2t> + 4g) + f | + 1- ) 



+ Q+|e 2 )sin(4i; + 4g)+^ + ^sin(5i; + %) + ^sin(6i. + 4^) + |^sin(7i; + 4 )? ) 



+ 0>(g), (6.33) 



where O(g) is a constant of integration. Here 



■§(*)W* 



(6.34) 



where 



q, = 3— 15^ + 



10? 



■i?; 



92 = 5i7 2 3 



:Vo 



qs 



35 



< 



(6.35) 



With use of (6.18) and (6.33) we find 



uL 7 __ 



M 4 +0(/ 2 ), 



(6.36) 



and 



*-!^(>-SD(>-'f) 

^ 3 " 64C? I 1 G 2 J' 



(6.37) 
(6.38) 
(6.39) 



through terms of order Jt 

To make Si purely short-periodic, one would have to choose 4>(g) in such a way that 
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27T.S, 



J 



S\dl would vanish. It turns out that <P(g) would then not vanish, but would have to be a 



long-periodic term, just cancelling a long-periodic term of order J\ arising from the rest of the 
expression for Si. The later calculation of long-periodic effects, however, will be accurate only 
through order J2. If we arbitrarily drop Qig), whose calculation would be extremely laborious, the 
net effect will be only to leave in the short-periodic terms of G, I, g, and h some long-periodic 
impurities of order /i, not affecting the accuracy of the calculation. 

The errors of the short-periodic terms will be of order /!, for two reasons. First, the short- 
periodic terms of the reference orbit were calculated only through terms of order J\\ second, the 
present calculation makes use of elliptic approximations, so that a variation of the form 0-4/* has an 
error of order J 2 in / and thus of order J\ in a^f. For this second reason the secular corrections 
produced by the perturbing potential will also have errors of order J% even though the first omitted 
term in the von Zeipel equation (6.25) is of order a\ or ]\. 

7. Short-Periodic Terms 



It is now straightforward but tedious to calculate the short-periodic terms. From (6.22), 
(6.31), (6.32), and (6.15), we find, through terms of order J\ 



dSi — r r > — ~ t 

—r — L> L — 0^k4, 



(7.00) 



^here 



^=-8^{( 3 - 15 ^ + f^)[(7) 5 - (1 - e2) " 7/2(1+ l e2 ) 



- 1 rfi) ("Vco»(2» + 2g) - ^ (1 - e 2 )- 7 ' 2 cos 2g 



+ -g-i7JFJ cos(4t; + %)j. 



(7.01) 



From (6.22) and (6.33), we find 



where 



dSi 



= G — G f — 0*4^4, 



(7.02) 



Ga = ° 2 [~ I e2{v " /} sin 2g ~ i cos {v " 2g) + ( 3e + t) cos {v + 2g) 

,(2i; + 2g) + (e + ^cos(3t; + 2g) + ^co8(4i; + 2g) + ^cos(5i; + ^)J 



■f ( 1 -r ~z e 2 ) cos i 



+ (?s 



| e 3 cos (v + 4g) + 1 e 2 cos (2v + 4g) + (2e + | e 3 ) cos (3v + 4g) 4- (1 + 1 e 2 ) cos (4t> + 4g) 



6e , 3e' 



1 



+ ( — + 777 cos(5i; + 4g) + - e 2 cos (6?; + 4g) + yj cos (7z; + 4g) 



5 , 10 ^_ v _ . . 0/ . 2 . w „ xw , w , 14 
Since the dependence of Sion H is only through the (?'s, the calculation of 



]■ 



(7.03) 



aSi _ 

h — h' — — -777 — 0-4/141 



(7.04) 



is simple. First calculate the derivatives dQj/dH. From (6.37) through (6.39) and the sufficiently 
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accurate relations H/G = cos /, H 2 /G 2 = 1 — tj 2 and G 2 = fip, we find 



3(?i = 15/r e 

a// 16 \p 

BQz_5(re\ 4 



dH 4\p 



(77J 2 — 4) cos/ 
(3 — 7t/ 2 )cos/ 



d<? 3 _35/r P y 2 . 



(7.05) 



From (7.04), (6.33), and (7.05) we then obtain 

-(^Ycos/{3(77?g-4) 
3 



16 Vp 
-h 4(3 - 7778) 



1+^e 2 1(»-0 + 



3e 3 \ 3e 2 e 3 

3e + — — Jsin v + — r- sin 2t> + 77- sin 3i; 
4 / 4 12 



e 2 (z; — /)cos 2# + — sin(i; — 2g) + ( — + -g- )sin(i; + 2g) 



3e , 3e ; 



1 , 3e 2 



e . e i 



3e 2 



(-+ — lsin(2i; + 2^) + (- + -lsin(3z; + 2^) + — sin(4^ + 2g) + 4 



+ 7^ 



3 -I 

- sin(5i; + 2#) 

V> : * V /p p*\ ( 1 V 2 \ 

I sin(i; + 4#) + -g- sin(2i; + \g) + f - + - jsin(3t> + 4#) + f - + — jsin(4i; + 4g) 

/3e 3e 3 \ e 2 e 3 

+ wo + 40 r' m{Sv + 4 ^ + s" sin(6x; + 4 ^ + 56 sin{lv + 4 ^ 



(7.06) 



To find /-/' and g-g\ we use (6.23), (6.33), and (6.38) through (6.40). We may write (6.33) as 



Si=o-4<2i(G, ff)/i(e, !;) + cr 4 <?2(C, #)£(*, *) + cr 4 <?a(G, //)/ 3 (e, t» f 



(7.07) 



wher 



/i(e,t;) = n+|e 2 Kt;-/)+^ (7.08) 



3 e' 5 /3e 3e 3 \ 

/2(c, *>) — 7 e 2 (*> ~~ ^) cos 2 # + ~q sin ^ ~~ 2 #* + \T + IT ) sin ^ + 2 #' 



^ (^ -r-^e^)sin(2^^ ^^T)^ ^^^^)sin(3^-r-2^)^^^ sin(4^ -r- 2^) + ^ sin(3^ -f- 2^) (7.09) 



/ 3 (e,t;)^^sin(t; + 4^)-h^sin(2t;+2g) + (|+^-j 



1 3 2 \ /3 Q 3\ 2 3 

-+-| : jsin(4t;+4g)+(^ + ^J S in(5 ! ; + 4g) + |sin(6i> + 4^ + |^sin(7i; + 4^. (7.10) 



Then 



dS 



dL °" 4 2<?'( de + dvde}dL 



fidv\ 

v de) 



de 



(7.11) 



i = \ 



dSj_ A dQif. + _ ft/)Yi6 + i£»is\i£. 



:<T <2^ + 



aC- 7 dL dG\dL 



(7.12) 
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From (6.37) through (6.39), with the same approximations used in obtaining (7.05), we find 

^ = l© 4(136 - 500r,§+385 ^ 

Next we need the elliptic approximations 

||=( / ita)-^e- 1 (l-e 2 ) (7.14) 

|J= (jZ^+j) s in v = (l -e 2 )" 1 (2 sin i;+|e sin 2t>) (7.16) 

To compute the dSi/dL and dSi/dG we need finally the six derivatives d/i/de and dfi/dv, i = l, 2, 3. 
From (7.08) through (7.10) we find 

d/i / 9e 2 \ 3 e 3 

■J- = 3e(t; — /) + I 3 + ~T") s i n f + 9 e sin 2v + — sin3i> (7.17) 

^=l+|e 2 + (3e + ^)cos?; + |e 2 cos2j; + ^cos3i; (7.18) 

U = | e(v - I) cos 2g+ ^ sin („ - 2g) + (| + ?f\ sin (« + 2g) 

3 /l v 2 \ 3*> 3e 2 

+|esin(2i; + 2gr) + (2 + T/ sin(3t ' + 2g ' )+ 8" sin(4j; + 2 ^ + ^ sin(5t;+2g) (7J9) 

2g+^cos (v-2g)+ (| e + ^) cos (» + 2*) + (l + |e 2 ) cos (2v+2tf) 

| e + yJ cos (3» + 2g) + ^- cos (4t> + 2g) + 1- cos (5w + 2#) (7.20) 



«/* 3 2 

— = -7e 2 cos 
dv 4 



^ = ^ sin (v + 4g) + j sin (2v + 4*)+(| + ^f) sin (3» + 4g) + 3 f sin (4» + 4*) 

+ ^+^) 8 in(5t;+4g)+|8in(6i; + 4g)+|^3in(7t;+4«') (7.21) 

|J = I cos (0 + 4tf) + ^ cos (2t> + 4*) + (y + ~) cos (3d + ig) + (l + y) cos < 4r + 4 ^ 

+ (y + t) cos (5r + 4g) + t cos (6t; + 4g) + 1 cos (7t; + 4g) ' (7-22) 

There then follow 

+ (~5 — tt) sin 3v + -r-sin4i;4- — sin 5» (7.23) 



(l-e 2 )(^ + U^=-|(14 + 5e 2 )sin2g+ 3 ^(l-e 2 )( l ;-/)cos2^ 

+ (t " If ) sin {v ~ 2g) + if sin (2u ~ 2g) + fi sin (3v " 2g) 
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\2 2 



19e 4 
16 



sin {v + 2g) + (| e - 1 e 3 ) sin (2d + 2g) + (| + e 2 - y^) sin (3v + 2g) 



m. 



^(^ + ^lsin(4i, + 2g)+l— + — lsin(5i; + 2g) 



'6e 2 , 3e 4 



"-«(f + f 



v de) 



Sp 3 p 4 

+ — sin(6i; + 2g)4 — sin(7i;4-2g) 



e + 



lie 3 
8 



l^ sin (1; - 4g) — jg sin 4g - M- 4- -^j sin (v 4- 4g) 
)sin(2i;44s)-(| + |e 2 4^^ 



/13 21 p 2 1V\ / 3p 3 \ 

+ (l§+ 15— if-J sin ( 5u + 4 #> + ( 2e + x) sin (6v + 4g) 

AV 2 9p 4 \ Sp 3 p 4 

+ (^- + §4) sin (7»+4*) + ^Mn<8B + 4g) + Jg sin (9v+4g)- 

From (6.34), (7.14), and (7.16) it follows that 



(7.24) 



(7.25) 



and then from (6.23), (7.11), and (7.23) through (7.26) that 



(i = l,2,3) (7.26) 



/ I — (T4I41 , 



(7.27) 



where 



(l-e 2 )* (r e y 



-)' 

PI 



hi = 2m (^) i 6 < 8 - 40 ^« + 35r >t) 



48(1 -e 2 )(i;-/)-f-(40+12e 2 -17e 4 ) sin?; 
e 



+ 4(20 - e 2 ) sin 2v + e(40 - e 2 ) sin 3v + 10e 2 sin \v + e 3 sin 5i; 



4^(7^-6) 



240(1 -e 2 )(i;- /) cos 2g4 5e 3 sin (Sv - 2g) + 50e 2 sin (2v- 2g) 

10, 



+ 15e(16 - 3e 2 ) sin (t; - 2g) - 20(14 4- 5e 2 ) sin 2g4 — (8 - 24e 2 - 19e 4 ) sin (v -f 2#) 

4- 240(1 - e 2 ) sin (2t> 4- 2g) + — (24 + 16e 2 - 5e 4 ) sin (3v + 2g) + 20(17 4- 2e 2 ) sin (4v + 2#) 



+ 3e(64 + e 2 ) sin (5v + 2g) + 50p 2 sin (6v + 2g) + 5e 3 sin (7» + 2g) 



+ t?3 



35e 3 sin (v - 4g) - 350e 2 sin 4g - 105e(8 + 5e 2 ) sin (v + 4g) 
70, 



- 140(8 + 1 le 2 ) sin (2v + 4g) - — (8 + 20e 2 4- 7e 4 ) sin (3v 4 4g) + 840(1 - e 2 ) sin (4i; 4- 4g) 

14 
4- — (104 4- 84e 2 - 13e 4 ) sin (5v + Ig) + 140(16 + 3e 2 ) sin (6v + 4g) 



4- 15e(88 4 3e 2 ) sin (7v 4- 4g) 4 350e 2 sin (8v 4- 4g) 4 35p 3 sin (9v + 4g) 
From (6.23), (7.12), (7.14), and (7.15), we find 

g-g'=~a-e 2 )-Hl-l')-cT 4 ^^fr 



(7.28) 



(7.29) 



699-888 0-63-2 
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Then, from (7.13) and (7.08) through (7.10), we find 

g-g' =04#«, (7.30) 

where 

g4l =-(l- e 2 )-* Ui ~ 3Y2 (^) 4 {2(136-500t ? ^ + 385^)[6(2 + 3e*)(v - 1) 

+ 9e(4 + e 2 ) sin v + 9e 2 sin 2i; + e 3 sin 3v] —2(12 - 82r?§ + 1 l7)3)[60e 2 (i; - /) cos 2g 

+ 10e 3 sin (v - 2g) + 30e(4 + e 2 ) sin (i; + 2g) + 20(2 + 3e 2 ) sin (2v + 2g) 

+ 10e(4 -f e 2 ) sin (3i> + 2g) + 15e 2 sin (4v + 2g) + 2e 3 sin (5» + 2g)] 

+ 7jg(l lryg - 4)[35e 3 sin (t; + 4g) + 105e 2 sin (2i> + %) + 35e(4 + e 2 ) sin (3v + 4g) 

4- 35(2 + 3e 2 ) sin (4v + 4g) + 21e(4 + 3e 2 ) sin (5i; + 4#) 

+ 35e 2 sin(6i; + 4g) + 5e 3 sm(7i; + 4g)]}- (7.31) 



This concludes the solution for the short-periodic terms arising from the residual fourth harmonic, 
with coefficient cr 4 —74+71- 

8. Long-Periodic Terms 

By (6.21) and (6.28) the Hamiltonian F* is 

F* = Fo(L',G f , H)+FUL\ G' ,H) + F lp (L' , G', //, #'), (8.00) 

short-periodic terms having been eliminated. We now try to find new canonical variables L", 
G", #", l\ g", and A", corresponding to a new Hamiltonian F**(L", G", H") + F**{L\ G", H"\ so that 

Fo(L\G\H)^Fi c (L\G\H)+F lp (L\G\H,g f )^Fr(L\G%H , )^Ff*(L\G\H f '). (8.01) 
If we can do so, then L", G", and //" will be constants of the motion and 

aF** „ aF** . aF** 
l= -^u ^=~^ h== —8lF' »«*> 

To find the necessary canonical transformation, we introduce the generating function 

S*=L , T+G"g'+H"h' +SW, G", //", g') (8.03) 

where Sf is to be of order cr 4 . Then 

j , d§_ rt, ifr = dS* _ j, dS* 

L ai' L L du dL" 

G= 17 =G+ ^ (8 ' 04) * = W> =g+ ]c; (8 - 05 > 

0_ rr/_ dS* _ „„ h»— — — W I d^l 
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Insertion of (8.04) and (8.05) into (8.01) leads to 

f.(l>,v+§ ,h) + f 1c (l>,g* + § >h )+f 1p (l>,g»+§ ,h,* + §, 

= F?* + Ff*, (8.06) 
whose Taylor expansion, with neglect of terms of order a\ or higher, is 

F4JL\ G",H)+T?%^T+FUL',G",H)+Fi P (L',G",H,g')=F$* + F?*. (8.07) 

oG ag 

Splitting (8.07) into zeroth order and first order terms yields, respectively, 

F$*=F (L',G'\H) (8.08) 

3§af + ^c + F lp = FT*. (8.09) 

Resolution of (8.09) into constant and long-periodic terms then shows that 

Ff*=F u (L\C\H) (8.10) 



gfo BSt 

dG" dg 1 



7,^=-F lp (L\G",H,g'). (8.11) 



With use of (4.17) and (6.13) and of double primes to denote quantities corresponding to G", (8.11) 
becomes 

2vr K - 1® f^ = ^^ e"H 1 - e" 2 )^ 2 r,'^ - \ iff) cos 2ff' . (8.12) 

By eqs (7.34) and (7.37) of an earlier paper [Vinti 1961 a], we have 

2n(v';- «#=-^r(5 cos 2 /-i)+oap 

J^f i-^ + otn (813 > 

4 L'*G"* \ G" 2 ) U2J ' 

On inserting (8.13) into (8.12) and replacing elements by Delaunay variables, we find 

AS* 5/x 2 r 2 o- 4 /, G" 2 \ / // 2 \ / H 2 \ / r ffV 1 

^ = "To!^ C ° S2 ^ » 14 > 

Integration then yields 

Then 

as* 5^0-4 / /7 2 W // 2 W // 2 \-> 

JU = -W^^A 1 'G^)V~ 7 G^){ l - 5 G^) * ln2g (al6) 

dG" 32y 2 G"4 \ G" 2 ) L' 2 \ G" 2 ) G" 4 \ L' 2 )\ G" 2 ) 

-80^(l- — )(l-5^i) 
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sin2g-' (8.17) 



dSf 5uV(7 4 ff / G" 2 \r H 2 I // 2 \-' H 4 / H 2 \~ 2 1 



We now express the long-periodic terms in the notation of (5.00) and (5.01). 

-G =-rG A g-g=-rg4 

Jl J2 

it in °"4 T 1, i„ °"4 r 

I — i — — 14 n — n — — AI4. 

J2 J2 



(8.19) 



If we then use (8.04), (8.05), (8.14), and (8.16) through (8.19) and place G" 2 = iip, L' 2 = fjLa, HjG" 
= cos 7, G" 2 /L' 2 =\ — e 2 ,g'~ unperturbed g as given by (4.18), we find to the required accuracy 



C 4 

7 4 = 



5r 2 ra 



*t e2 ( l ~ e*)- 3l2 7)&l - 7 cos 2 7) (1 - 5 cos 2 1)- 1 cos 2g 
16 



16 U/ 



(1 - e 2 )~* r)l(l - 7 cos 2 /) (1 - 5 cos 2 7)- 1 sin 2g 



(8.20) 
(8.21) 



g4 = -^2 (~) 2 [2 + e 2 - 3(2 4- 3e 2 ) cos 2 / - 8(2 + 5e 2 ) cos 4 7(1-5 cos 2 7)" 1 

- 80e 2 cos 6 7(1 - 5 cos 2 7)~ 2 ] sin 2g (8.22) 

h 4 = - ^ (-Ye 2 cos 7[3 + 16 cos 2 7(1-5 cos 2 7)- 1 + 40 cos 4 7(1 - 5 cos 2 1)' 2 } sin 2g. (8.23) 
16 \pj 

9. Secular Effects 



We now have to use (8.02), (8.08), (8.10), and (4.17) to obtain /", g", and h" as linear functions 
of the time. From (4.17) and (8.08) we obtain 



dF$* = o „ sf$* 



a/'?* 



BL , —1 dC „ ■■2tt{v'1-v'£ -^-=27r(^sgna 3 -^. (9.00) 



From (8.10) and (6.12), with use of e' 2 = 1 -G" 2 /L' 2 , L' 2 = fxa\ and ^ 2 - l-77 2 /G" 2 , we have with 
sufficient accuracy 



Then 



/T**: 



3/x 6 rfo- 4 
128L' 3 G" 7 



/ ' 2 / 1 f~"' 2 /^'m j * 



3Ff*_ 45/^4 / G» 2 W 77 2 77 4 

dL' 1287.' 4 G" 7 V 7>' 2 / V G" 2 G" 4 

dF** 15/LtVo-4 



dG" 128L' 3 G" 8 
aFf*_15/x 6 r^o- 4 



G" 2 \H 4 



( G" 2 \ I G" 2 \ H 2 l G" 2 \ 

3^7-3— j + 18(^7 -^-15j —+7(^55-27 -^j G , 



677 32L' 3 G" 



>-3^)(3-7^)- 



(9.01) 

(9.02) 
(9.03) 
(9.04) 



Since our constant orbital elements a, e, and 7, and thus the frequencies V\ and i>2, correspond 
to 7/, G", and 77, we may drop the double primes from the i>'s in (9.00). Also, to our specified ac- 
curacy, we may put e 2 = 1 -G" 2 /L' 2 , L' 2 = ^a, and r) 2 = l-H 2 IG" 2 in (9.02) through (9.04). Then, 
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with the unperturbed values /, g, and h given by (4.18), we find 

/" = / + 0-4/42 g" = g + (T4g42 h" = h + 0-4^42, 

where 

/42 = -^(-)%l-e 2 )^(3-30cos 2 /+35cos 4 /K 

g42 = ~~ ik (") ^ ^ 3(4 + 3e2) ~ 18(8 + 7e2) cos2 1 + 7(28 + 27e2) cos4 /] / 



h\2 — ~ 



15 /ft 

32 



£)■ 



(2 + 3e 2 )cos/(3-7cos 2 /)ft 



(9.05) 

(9.06) 

(9.07) 
(9.08) 



We postpone summarizing the results for the residual fourth harmonic until the complete algo- 
rithm in section 11. 



10. Effects of the Third Zonal Harmonic 
By (2.01) the contribution of the third zonal harmonic to the potential is 

AV= t irlr- 4 J: i P,(s\n6). 
Corresponding to (3.02), this leads to 

AFi=-fifJr-V»P«(rinfl) 
and thus to 



Aft--*?/. 

a 4 



;Vo + -^Vo(z) sin (t> + #)- 5 *??(-) sin(3i; + 3g) 



8 



8 



corresponding to (6.06). Following the notation of section 6, we then find 



which is purely long-periodic, so that 

AF lm = AF lp (10.04) AFi c = 0. 

Since AFi c = 0, there are no secular effects. 
The short-periodic part of AFi is then 



AFn = AF 4 -AFi 



fir* e 



'2 Vo *T V v\\r) sin ^ + ^ 



(10.00) 



(10.01) 



(10.02) 



AFi w = (2tt)- j j^ AFK// = -^y3e(l-e 2 )- 5 / 2 (-|T ? o + y^sin^, (10.03) 



(10.05) 



- e(l - e 2 )~ 5/2 sin gl - 1 rft ("\ sin (3v + 3g)l • (10.06) 
On following the same procedure as in section 6, we find 



dAS l 



= -(2tti/i)- 1 AFu 



(10.07) 
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where 



V>rg 



ASi--(2wi)- 1 /AFnd/ 
; 73(1 - e 2 )- 5 / 2 T/ {3(l - 5 cos 2 1) [ 12e(i; - /) sin g + 3e 2 cos (v - g) 



(10.08) 



— 6(2 -he 2 ) cos (v+g) — 6e cos (2t> + g)— e 2 cos (3v + g)] -h T7§[15e 2 cos (t> + 3g) + 30e cos (2i;-h3g) 

+ 10(2 -I- e 2 ) cos (3v + 3g) + 15e cos (4v + 3g) + 3e 2 cos (5i; H- 3g)] }. (10.09) 



Before taking derivatives of (10.09) with respect to V and G' one must replace e 2 by 
1 — G' 2 /Z/ 2 , 7] 2 by 1 — H 2 IG' 2 , and a by L' 2 I/jl. The short-periodic effects are then given by 



aASi 

a/ 



=L-L'=JsL 3 



a as, 



___ G _ G ' =i /3 G 3 (10.10) -^■ = g - g '=J 3g3 



^here 



Ls = - 



nr? 



_aASi_ 



-2 T ?o + ~8~ T >o 



H-H' = 



dASi 
dH 



— h — h'—Jzhz, 



— ) sin^ + g) — e(l — e 2 ) _5/2 sing 



-pi r sin(3i;-h3g) 
G 3 = -na 3 (l-e 2 )^(— j K""^ 170 ^'^^) e ( t; ""^) cos ^+ ^sin(t;-g) 



(10.11) 



(10.12) 



1 \ 1 e 2 

+■ ( 1 4- - e 2 1 sin (v + g) + - e sin (2t> + g) + — sin (3*; + g) 



15 fe 2 1 

— 178 — sin (v + 3g) +- e sin (2v + 3g) 



+ 



|(l+|e 2 Vsin(3t; + 3g)+|sin(4t;-h3g) + ^s 



20' 



(10.13) 



/3 = g(l-e 2 )*g) 



( 1 — e 2 ) (v — I) sin g -f I 1 + — e 2 J cos g 



e 1 <? 3 

-j- — (14 — 3e 2 ) cos {v — g) + -e 2 cos (2i? — g) + — cos (3i; — g) — e(l — e 2 ) cos (v + g) 
o z lo 

3 e 1 e 3 

— - cos (2v + g) - 7^7 (34 — e 2 ) cos (3v + g) - - e 2 cos (4i> + g) — — cos (Si; + g) 



■57|g 



1 e 3 e /l \ 

-e 2 cos 3g — — cos (v — 3g) — -(6 + 5e 2 )cos (^ + 3g)— (--he 2 1 cos(2i;H-3g) 



_±„2 

2 



+ 1 (1 - e 2 ) cos (3v + 3g) + 7 (5 + e 2 ) cos (4>v + 3g) -f ^r (54 + e 2 ) cos (5v + 3g) 
o 4 40 

1 e 3 

-f - e 2 cos (6f + 3g) -h 7-7 cos (7i> + 3g) 
z lo 

g 3 = -(l-e 2 )-^/ 3 -|^) 3 {r ? o 1 (4-39^ + 40^) e(v- /) sing + ^ cos (v-g) 



(10.14) 



1 \ 1 e 2 1 

1 + - e 2 1 cos (t> + g) - - e cos (2v + g) — — cos (3v + g) + 5T7o(8-n 2 - 3) 



■ cos (v + 3g) 
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+ 2 e cos (2v + 3g)+ x (l +1 c 2 ) cos (3t; + 3g) + 1 cos (4v 4- 3g) + ~ cos (5t; 4- Sg) 
^3 = (~) l^cot/l 1 — — sin 2 / j e(v — /) sing- 4- — cos (v — g) — f 1 -\-~e 2 \ cos (v + g) 



(10.15) 



1 e 2 

— - e cos (2v 4- 2g) - — cos (3v 4- 2g) 



— — sin 11 
lo 



— cos (v + 3#) 



1 1 / 1 \ e e 2 

4- - e cos (2v 4- 3#) 4- - ( 1 + - e 2 1 cos (3t> + 3g) 4- - cos (4v 4- 3g) H- — cos (5v 4- 3#) 



(10.16) 



Continuing on to the long-periodic effects and following the procedures of section 8, we find 

(AFi)** = 0, (10.17) 



so that the third harmonic gives no secular changes, and 



(10.18) 



Now 



, 3 ^r'je" r n (,,H?_ 

lp or tzfnsJ&HQ M 3 ^tt 2 



(10.19) 



by (10.03) and (10.04). With use of (10.18) and (10.19) and of (8.13) for 2tt(v'[-v'^, we find 



so that 



dASf _l fireJ-A ,, „ . , 






(10.20) 



(10.21) 



It follows that 



dl' ~ L L ~ u 

dAS* 



dAS* 



dh 



t-=H'-H"=0 






-l'-l"=hklh 



dAS* 



G'-G" = J,GM 2 (10.22) — 3^r-=«r'-g" 



dAS* 
Ojjr — h' — h" = JMJi, , 



(10.23) 



where 



G,\ = - r,na(\ — e'-) * erasing 
/:, = o p sin/eosg 



g3 z 



1 r e le cos 2 / sin / 



2 p\ sin/ e 

hz = — - — e cot / cos g. 
2p 



cos g 



(10.24) 
(10.25) 

(10.26) 

(10.27) 
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A summary of the results for the third harmonic will appear in the next section, which gives the 
complete algorithm. 

1 1 . The Complete Algorithm 

We shall summarize results by writing out the complete algorithm for the calculation of the 
motion in the potential field (1.00), as modified by the third harmonic and the residual fourth 
harmonic. This will involve a repetition, with some changes, of about two pages of an earlier 
paper [Vinti 1961 a], but it is highly desirable to assemble the whole solution in one place. 

Let the planetary constants p, r e , J \, J 3, and J 4 be given, along with the constant orbital ele- 
ments a, e, /, /o, go, and j8 3 . To calculate the 



compute 



Unperturbed Reference Orbit 

c = r e j£ 170 = sin/ p = a{\ — e 2 ) 

D = (ap- c 2 ) (ap - c 2 T)l) + 4a 2 c 2 7? 2 ) D' = 4a 2 c 2 (l - tj 2 ) + D 

A=- 2ac 2 D~ 1 (l - 7) 2 ) (ap - c 2 r) 2 ) B = c 2 r) 2 D- l D' 

bi=~~A b 2 = B* -2a 1 = p(a + b 1 )- 1 

~2 = a po = - c 2 (1-t? 2 ) + apD-W a 2 = (-2ai)* (a*p )* >0 

a3=a T"W cos/ % > ~^v k=cp q=r}oV > 

A^d-e^p £ (b2lp) n Pn(bi/b 2 )Rn-2[(l-e 2 )*] 

n=2 

A 2 = (\-e*)*p-' J (6 8 /p) B P»(fci/62)/?»[(l-e 2 )'*], 

n=0 

where P n (x) is the Legendre polynomial of degree n and where R n (x) = x n P n (x~ l ), always a poly- 
nomial of degree [n/2] in x 2 . 

A 3 = (l-e 2 )*p-* J D m /? m+2 [(l-e 2 )*], 

m=0 

where 

D 2i = V ( - l) i - n (clp) 2i - 2n (b2lp) 2n P2n(b l lb 2 ) 
n=0 

D 2M = 2 (-iy- n (c/pf i - 2n (b 2 lp) 2n+, P2n + i(bilb 2 ) 

n=0 

B l = 27T-'q- 2 [K(q) -E(q)] £ 2 = 27T- 1 K(< ? ), 

where K(q) and E(q) are the complete elliptic integrals of the first and second kinds, respectively. 
It may be convenient to have the series 
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Continue with 



R,= l— (1 _„-*>-* _ V V T)" 2 " 1 where v _ ( 2/ra )! ys' ^n)!?) 2 ," 

° 3 L u "* > 2, Y m1 ?2 ' wne re y™ 2 2m (m!) 2 ^ 2 2 "(n!) 2 

^„=|(l- e 2 )V- 3 e(- 26,6^ + 6^) ^,2 = ^(1 -e 2 )V" 3 ^ 2 



^2i = (l-e 2 ) i4 p- 1 e 
^22 = (l-e 2 )^p- 1 



6,p- 1 + (36f-6|)p- 2 -|6 1 6l(l + 



iy+h 



(4 + 3e 2 )p- 4 



| (36 2 - 6 2 )p- 2 - 1 e 2 Mfp- 3 + ^ 6 4 (6e 2 + e«)p-< 



^23 = (1 - e 2 ) V 1 1 ( - M|p- 3 + btp-*) 



A i4 = ^(l-e 2 )*p-*bie* 



A 3l =(l-e 2 )*p- 3 e 
A 32 = (l-e i )*p- 3 
A xx = (l-e 2 )*p- 3 e 3 



2 + 4,p-'3 + 



3e ; 



)-„-i\ 



4 / 

4 2 



6 2 + c 2 )(4 + 3e 2 )| 



f + f^p-^-p-^ + fe 2 )^ 



+ c 5 



^34=-^(l-e 2 )V- 5 e 4 (|^ + c 2 ) 

2tt^ 2 = (<*! - a§* tfo x A&i\a + 6, + y4i + c 2 7ilA 2 B x B^)-\ 

The uniformising variables £, t>, and i// are then given by E = M s -\-E p , v = M s + E p , and i// = i// s 
+ i// p . If £ is the time, their secular parts M s and i//. s are given exactly by 

M s = /o + 2lTPit l//,- = /o + #0 + 2771/2^. 

Let the periodic parts be split as follows: E p = Eo + E\+E2, v p = Vo + v\ + V2, and i// p = i//o + */>i 
+ 1//2, where, e.g., E contains terms of order /§, J 2 , and /§, Ei contains terms of order J 2 and 7|, 
and ^2 contains terms of order J\ only. 

Then Eq is given by the Kepler equation 

M 8 + £ - e' sin (M, + £ ) = M„ 

where e' = ae(a H- 6i) _1 < e. The term v o is then given by placing t> = M« + v o and £ = M. s + 1; o in the 
anomaly connections 

cos v = (cos £ — e) (1 — e cos E)' 1 sin f = (1 — e 2 )* (1 — e cos £')" 1 sin £ 

or equivalent relations. (Note that e here is the original e and not the e' in the Kepler equation.) 
Then 



The term E\ is now given by 



«/*> = (- 2a,)"* (oi-al^ijo^iflrH*. 



1 , r 



£ 1 = [l-e'cos(M s + £o)]-W 1 -|eTl-e'cos(M s + ^o)]- 3 M?sin(M s + £ , ), 
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where 



M^ia + biY 1 



-(Ai + c^^BiB^vo + ^ (- 2aO* (af - <*§)"* T7§ sin (2<fc 4 2i// ) 



The term vi is then given by placing v = M 8 + v + vi and £ = M s 4 £ 4 £1 in the anomaly connec- 
tions. Then 

i/n = (- 2ai)-^ (a|- a§) V4 T7o 1 B2 1 [^2Vi +^21 sin (M, + 1> ) + ^22 sin (2M S 4- 2i> )] 

-h^- 1 sin(2i// s + 2i// ). 



Finally 



/here 



^^[l-e'cos^ + ^o + ^i)]- 1 ^, 



M 2 =-(a + 6i)" 1 L,i;i+^nsin(M,-ht;o) + ^i2sin(2M. s + 2i;o) 



+ c 2 (- 2«i)* (a| - «§)"* ijg j £h|*i - |^i cos (2i// s 4 2i// ) - ^ sin (2i// s + 2i// ) + ~fe sin (4t// s 4- 4i// ) 

Then t* is found by placing v = M 8 4 v 4 vi 4 v 2 and E = M 8 + Eo + Ei + E 2 in the anomaly connec- 
tions and 

1//2 = (- 2ai)-* (af— a 2 .)* 'no l Bi l [A%V2 +A21V1 cos (M« 4 v ) 4 2^ 2 2i>i cos (2M, 4 2t; ) 



4 ^23 sin (3M, 4 3i; ) 4 ^ 24 sin (4M. S . 4 4v )] + 7^ 



i// lC os(2i// s 4 2i// ) 



+ ^ sin (2i//, 4 2i/*>) - |£ sin (4i// s + 4i// ) j 



Unperturbed Reference Orbit 
The spheroidal coordinates p and rj are then given by 

p = a{\ — e cos£) = (l 4e cos v)~ l p V = r )o sin 1//, 

where E = M s -\- Eo~\- Ei+ E 2 and v = M s 4 ^0 4- vi + ^2. To obtain the unperturbed right ascension <£, 
first calculate an angle x> equaling 1// whenever xfj is a multiple of 7r/2 and satisfying 

cos x = (1 ~~ Vo sin 2 i//) - ^ cos 1// sin x = (1 — Vo)* (1 ~~ '"o sin 2 ^) - ^ sin 1// 

Then 

(/> = ^4a 3 («i-a 2 3)-^T7o 



(l-T ? 2 )-Ml-^ 2 - 2 )-^X + ^3t// + ^r ? 2 T ?2 - 4 sin2i// 



] 



- c 2 a 3 (- 2«i)-^ \A 3 v 4 2 ^ 3 * sin 



The Zonal Harmonic Perturbations 

Within the accuracy of calculation of the perturbations, we may use either the approximate or 
the accurate formulas listed below. (Use of the approximate formulas may involve more work than 
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use of the accurate ones, because the latter will already be known from the solution of the un- 
perturbed problem.) Compute 





Approximate formula 


Accurate formula 


n 
I 


/o + nt 


Mg = l + 2irvit 


g 


go + §(^)V 2 (5cos 2 /-l)* 


gQ + 27T{v 2 — V\)t 


E 

V 


M g + E 

M s + Vo 


Ms + Eo + Ei+E-, 
M s + v + vi + v-> 



Also, r = p, within the accuracy of the calculation. Then compute 

Third Harmonic Short-Periodic Quantities 



t 3 =- 



(-|i7o + ^i72) {(f) 4 sin (f+^)-e(l-e 2 )- 5/2 3ing}-|r,g(;) 4 sin (3t;+3g) 



Gz = — na 2 ( 1 



-^ egr(-M« 



e(v — l) cosg-\~ — sin (v—g) 
4 



1 \ 1 e 2 

4( i + o 6 / sin ^+#)4-esin(2i;4g) + — sin 



15 



jsin(i; + 3^) + |esin(2i; + 3^) + |(l+|e 2 



(3i> + *)] 

I sin (3v 4 3#) 



4- 7 sin (4v + 3g) + — sin (5v 4 Sg) 



]}• 



(Note that v — l — vo, within the accuracy of the calculation.) 



i»=f e a-er 



( & W5i)5-4) (l-e 2 )(i;-/)sin^+(l + |e 2 )cos^4|(14-3e 2 )cos(i;-^ 



1 e 3 3 e 

4-e 2 cos (2i; — g) 4 77 cos (3z; — g) — e(l — e 2 ) cos (v + g) — -cos (2i;4#) — 77(34 — e 2 ) cos (3i;-r-g) 
z lo z z4 



1 e 3 

■ - e 2 cos (4f 4- g) — 77 cos (5v 4- g) 
z lo 



+£<!-«■>« 



fcM-s' 



2 cos 3^ — 77 cos (v — 3g) 



- | (6 + 5e 2 ) cos (v + 3g) - ( 1 4- e 2 j cos (2?;-f3g) + |(l -e 2 ) cos (3i; + 3g)4|(5 + e 2 ) cos (4i;4 3g) 

4 ^r (54 4 e 2 ) cos (5» 4 3g) 4 1 e 2 cos (6v 4- 3g) 4 ^ cos (7v + 3g) 



#3 



= - (1 - e 2 )~* / 3 -| (*f) V<4 - 39t, 2 + 40V ) 



e(f — /) sin g 4- 7- cos (t> — g) 
4 



5 (r e ^ 



cos (t; +g) — 2 e cos (2t; 4 g) — — cos (3?; 4 #) 



^(8^-3) 



e 2 1 1 / 1 \ 

— cos (v 4 3g) 4 - e cos (2?; 4 3g) 4- - I 1 4 - e 2 1 cos (Sv 4- 3g) 



4- 7 cos (4t; 4- 3g) 4- — cos (5v 4- 3g) 
4 zU 
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fe -(?) , {l c »"( 1 -f"0 



e 2 / 1 

e(i; — /) sin g + — cos (t> — g) — ( 1 + - e 2 ) cos (v 4- g) 



1 



- - e cos (2v + g) - — cos (Sv 4- g) 



12 



15 . 9 , 

- — sin 2/ 
lo 



e 2 1 

— cos (v 4- 3g) 4- — e cos (2t> + 3#) 

4 2 



+ 



t(>+B 



3 \* ' 2 ~ / w " v ~" ' ~ e/ ' 4 
Third Harmonic Long-Periodic Quantities 



cos (3v + 3g) + - cos (4>v 4- 3g) 4- — cos (5i; -f 3g) [• 



Gz = -f e na{\ — e 2 ) ^ erasing 



/a = ^ (1 - e 2 )** Vo cos g 
r e (e cos 2 / sin A 

g3 — 7T : T COSg 

2p \ sin / e J 

hz=—-£-e cot / cos g 



,- "l(a) 



(/ua)» 



Fourth Harmonic Short-Periodic Quantities 
3-15T,J+^^)[^) 5 -(l-eV /2 (l+|e*) 



+ 5^(3-|^)[(f) 5 



cos (2v + 2s)-^e 2 (l -e 2 )" 7/2 cos 2g 



+ ~8 V °\r) cos(4l; + 4 ^ 



^=-|(f) 4 (MP)^s{(3-|r, ? ) 



3 e 3 

— - e 2 (f — /) sin 2g — — cos (v — 2g) 
2 4 



3e + ^)cos(i; + 2#)+(]+|e 2 )cos(2i; + 2g)+(e + - 



+^cos(4t;+2^)+^cos(5t;+2 



] + I^o 



1 3e 2 

- e 3 cos (v 4- 4g) 4- — cos (2t> 4- 4g) 



2c4-|e 3 ) cos (3v4- 4#)4- (l+f* 2 ) cos (4z; + 4#) 4- (y + ^) cos (5z; + 4#) 



' 41 1 P ; 2Q48 



1 e 3 

4- - e 2 cos (6i; 4- 4>g) 4- — cos (7v 4- 4#) 

2 14 



6(8-40^ + 35^) 



48(1 -e 2 ) (*;-/) 



4- - (40 4- 12e 2 - He 4 ) sin v 4- 4(20 - e 2 ) sin 2v 4- e(40 - e 2 ) sin 3i; 4- 10e 2 sin \v 
e 



-he 3 sin 5i> 



- 4r? 2 ) (77) 2 ) - 6) |"240( 1 - e 2 ) (v - 1) cos 2g 4- 5e 3 sin (3i; - 2g) 



4- 50e 2 sin (2v~ 2g) 4- 15e(16 - 3e 2 ) sin (i; - 2g) - 20(14 4- 5e 2 ) sin 2g 

+ — (8 - 24e 2 - 19e 4 ) sin (v 4- 2#) 4- 240(1 - e 2 ) sin (2v 4- 2g) 4- — (24 4- 16e 2 - 5e 4 ) sin (3t> 4- 2g) 
e e 



4- 20(17 4- 2c 2 ) sin (4t> 4- 2g) 4- 3e(64 4- e 2 ) sin (5v 4- 2g) 4- 50e 2 sin (6v 4- 2g) 4- 5e 3 sin (7t> 4- 2g) 
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vi 



35e ;i sin (v - 4#) - 350e 2 sin ^g - 105e(8 + 5e 2 ) sin (v + 4g) - 140(8 + 1 le 2 ) sin (2*; + 4g) 



70 14 

- — (8 + 20e 2 + 7e 4 ) sin (3t> + 4#) + 840(1 - e 2 ) sin (4v + 4g) + — (104 + 84e 2 - 13e 4 ) sin (5i> + 4g) 
e e 



140(16 + 3e 2 ) sin (6v + 4g) + 15e(88 4- 3e 2 ) sin (7» + 4g) + 350e 2 sin (8v + 4#) 



+ 35e 3 sin(9t; + 4g') 



», = - (1 - e 2 )"* / 4 i - iL ^ {2(136 - 500t? 2 + 385t$ [6(2 + 3e 2 ) (v - 1) 



+ 9e(4 + e 2 ) sin v + 9e 2 sin 2v + e 3 sin 3i;] - 2(12 - 82r? 2 + 77t? 4 ) [60e 2 (i; - /) cos 2g 
+ 10e 3 sin (v - 2g) + 30e(4 + e 2 ) sin (v + 2g) + 20(2 + 3e 2 ) sin (2v + 2g) + 10e(4 + e 2 ) sin (3v + 2g) 
+ 15e 2 sin (4v + 2g) + 2e 3 sin (5v + 2$]+ 7/8(1 l^g - 4) [35e 3 sin (v + 4#) + 105e 2 sin (2w + 4g) 
+ 35e(4 + e 2 ) sin (3v + 4g) + 35(2 + 3e 2 ) sin (4v + 4g) + 21e(4 + 3e 2 ) sin (5t> + 4g) 

+ 35e 2 sin (62; + 4g) + 5c 3 sin (7v + 4g)] 1 
1 + — e 2 I (0 — /) + ( 3e + 7 e :} I sin t> -h - e 2 sin 2t> 



A 4 i=-^(-) 4 cosi 



-f — sin 3*; 



+ 4(3-7r? 2 ) 



-^ 



3e , 3e :} 



+ 



e 2 (i; — /) cos 2g + — sin (t; — 2g) + I — + - 



3e 2 



J sin (v -h 2g) 
(|+^)sin(2i;+2g)-h(|e+^sin(3t;+2g) 

7i7 , ■ \2 V " 8/™ v ~" ' *°' 

3e 3 \ . 
4oJ Sm 

Fourth Harmonic Long-Periodic Quantities 



+ — sin (4v + 2ff) + ^ sin (5v + 2g) 
?3 sin (0 + ig) + ^ sin (20 + 4#)+ ( i e + Q sin (30 + ig) + (7 + 2g-| sin (4e + 4#) 

+ fg-f— Jsin(5t;+4^)+|-sin(6t;-h4^) + ^ si 



8 



5r 2 n 
~16~ 



G 4 = - : ^e 2 (l-e 2 )- 3 / 2 T ? 2 (l-7cos 2 /)(l-5cos 2 /)- 1 cos2^ 



-£©"<•■ 



A to' 



- e 2 )~* <n 2 (l - 7 cos 2 /) (1 - 5 cos 2 7)- 1 sin 2g 



g4 = ~ 32 l~i ^ + e " " 3(2 + 3e2) C ° S2 1 ~ 8(2 + 5e2) C ° s4 /(1 " 5 C ° s2 /)_1 

- 80e 2 cos 6 7(1 - 5 cos 2 1)~ 2 ] sin 2g 
h 4 = ~ j7 (-Ye 2 cos 7[3 + 16 cos 2 7(1 - 5 cos 2 7)" 1 + 40 cos 4 7(1-5 cos 2 1)~ 2 ] sin 2g. 

Fourth Harmonic Secular Quantities 
/ 4 2 = -^(-Vn(l-e 2 )^ (3 -30 cos 2 7 +35 cos 4 1)t 

g42= ~vk (~) 4 "t 3(4 + 3e2) " 18(8 + 7e2) cos2 1 + 7(28 + 27e2) cos4 * 

^42 = - 1| (-Yn(2 + 3e 2 ) cos 7(3 - 7 cos 2 /)*. 
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With cr 4 — 74+71, next compute the variations in the Delaunay variables 

SL = J 3 L 3 + (r 4 L 4 8H = 
8G = J3G3 + y G-s + 0-4G4 + -7- G 4 

8/ = ^3/3 + y ^3 + o- 4 (/4i + 142) + ~rU 

J2 J2 

8g = j3g3 + -rg3 + (T4(g41+g42) + -rg4 
J2 J2 

8h = 73^3 + 7- ^3 + 0-4(^41 + A42) + 7- A4 

J«2 J2 

The variations of the elements a, e, and 7/0 are then 

8a = — 8L Se=— 8L - (ae)-^-)* 8G 8t) = ^—^8G 
an jme \fi/ fJLrjop 

The variations in the uniformising variables £", v, and 1// are 

8£ = (p/a) (8/ ■+■ sin £8e) 

8i; = (l-e 2 )* (p/a) [8£ + (l-e 2 )" 1 sin £8e] 

8\fj = 8v-\- 8g 
The variations in the spheroidal coordinates p, 77, and </> are then 

8p — (pla)8a — a cos £8e + ae sin E8E 

8r) = {rj/r)o)8r) + 770 cos 1//81// 

80 = 8A + (l-7yg s i n 2^)-i cos /[8^-|( /[jl p)-vi sin2i//8G] 

The final rectangular coordinates X + 8X, Y-\-8Y, and Z -\- 8Z are then given by 

X+^+i(Y+8Y)=[(j> + 8p) 2 + <?]*[l-(r) + &r)y]* expi(4 + 8d>) 
8Z = p8i7 + T78p. 

12. Discussion of Results 

Since the chosen accuracy of the quantities appearing in the present perturbation, with factor 
of order 7i> is only that of an elliptic approximation, the variations in the Delaunay variables may 
be checked against those found by Brouwer [1959], by Garfinkel [1959], or by Kozai [1962]. 
Comparison of the above results with those of Brouwer or of Garfinkel shows that the long-periodic 
effects of the third and fourth harmonics and the secular effects of the fourth agree with Brouwer's 
or Garfinkel's, provided that one replaces their J 4 by 74+71- Comparison with Kozai shows that 
the short-periodic effects agree with his. 1 Similarly one can read out of Kozai's paper the long- 
periodic effects of 75, J7, and J9 and the long-periodic and secular effects of J% and 7s; in so doing 
one ought in principle to replace his 76 by 76 — 7f an d his 7s by J% + 7|» Dut tn -is w °uld De going be- 
yond the accuracy of the present calculation. Since the author's orbital elements differ from 
Kozai's by terms of order 72, the agreement with Kozai holds only through terms of order 72 for 
long-periodic effects and through terms of order }\ for short-periodic effects. 



1 On his page 451, however, in the first line for AG, the expression cos 2g should be cos g. 
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To compare accuracies, we construct the following table, noting that the author's reference 
orbit accounts for about 99.5 percent of the deviation of the earth's potential from spherical sym- 
metry. Thus my perturbation potential is only about 0.5 percent of Kozai's. 



Effects of 99.5 percent of deviation from sphericity 



Kozai 
Author 



Secular accuracy 

Through J\ 
Exact 



Short-periodic accuracy 
Through J\ 
Through J\ 



Long-periodic accuracy 

Terms do not exist 
Terms do not exist 



Effects of remaining 0.5 percent of deviation from spherity 



Kozai 
Author 



Secular accuracy 

Through J'.\ 
Through ]\ 



Short-periodic accuracy 
Through J\ 
Through J\ 



Long-periodic accuracy 
Through J\ 
Through J> 



Thus the advantages of the author's treatment are the exact solution for the secular effects 
arising from 99.5 percent of the aspherical deviation and the much shorter algorithm. The prin- 
cipal advantage of Kozai's treatment, arising in connection with the remaining 0.5 percent of the 
aspherical deviation, is his more accurate solution for the long-periodic terms. 

The present solution, like all previous perturbation theories, gives rise to the resonance 
denominator 1-5 cos 2 / in some of the long-periodic terms. These terms are thus not reliable 
if one considers inclinations / sufficiently close to 63.4° or 116.6°. For such inclinations one could 
improve the accuracy by boldly dropping the long-periodic terms with coefficient J 4 -\-J\ or, better, 
by superposing on the present treatment Izsak's [1962] solution of the problem of the critical 
inclination. 

The element e occurs in the denominators of 8e, Z 3 , g A , lz, #3, Uu and g 4 \ and thus also in the 
denominators of 8/, 8#, 8E, 8v, and 8t//. No corresponding trouble occurs in the coordinates, 
however. To test this point, reject all terms except those containing e~ l . One then finds / 3 = 

— £3, l3 = — g3, Ui= — g4i, I»3 = Gs, L 4 = G 4 , and G 3 = G4 = 0. Then 8e = 0(e°), 8E = 8v = 8l, and 
81// = Sv + 8g = 81 -h 8g = 0. Trouble could occur in 8r) or in 8(f) only through 8t//, so that 8r) and 8(f) 
do not become infinite for e = 0. Similarly, trouble could occur in 8p through the term 

— a cos E8e, which, however, does not become infinite, or through the term ae sin E8E. But 8E = 81 
=0(e~ 1 ), so that this term also remains finite. 

The quantity r) = sin / occurs in the denominators of 8r) , #3, h 3 , gs, and A3, and thus also in the 
denominators of 8g, 8h, St//, and 8(f). Again, however, the coordinates remain finite when sin 7 = 0. 
To test this point, reject all terms except those containing -n^ 1 = esc /. Since G3 and G 3 are of order 
170 and Ga and G 4 of order 77^, it follows that 8G is also of order 170 and 8170 = (1 — y)l)8G/(r)otJLp) remains 
finite. Thus the term sin 1//8770 in 817 remains finite. We also have #3 = ~ h :i = 0(r)^ 1 ), g :i = \ er e p~ l 
X cos 2 I esc / cos g,h 3 = — 2 er e p~ l cos / esc I cos g, and 8^j = 8g= 0{t\q 1 ). Thus the term 170 cos 1//81// 
in 817 also remains finite, so that 81) remains finite. Finally, 8(f) = 8h-\- cos I8g = Js(gs-\-h3)-\-JsJ~ l (hs 
-\-g3 cos I) = ^j3J~> l er e p~ l cos g esc / (cos 3 / — cos l) = — ?j3J2 i er e p~ l sin/ cos / cos g=0. Thus 
8(f) remains finite. There is no trouble with 8p. 



221 



13. References 

Brouwer, D. 1959, Astron J. 64, 378-397. 

Garfinkel, B. 1959, Astron. J. 64, 353-367. 

Garfinkel, B. 1960, Notes of the Summer Institute in Dynamical Astronomy at Yale University, July 1960, page 281. 

Izsak, I. 1962, Smithsonian Astrophysical Observatory, Special Report No. 60. 

Kaula, W. 1962, Advances in Geophysics, Vol. IX, edited by H. E. Landsberg and J. Van Mieghem, p. 225 (Academic 

Press). 
King-Hele, D. G., C. E. Cook, and J. M. Rees, 1963, Nature 197, 785. 
Langer, R. M. and J. P. Vinti, (1963) J. Res. NBS 67C (Eng. and Instr.) No. 3, 247-249. 
Vinti, J. P. 1959a, Phys. Rev. Letters 3, 8; 1959b, J. Res. NBS 63B (Math, and Math. Phys.) No. 2, 105-116; 1961a, J. Res. 

NBS 65B (Math, and Math. Phys.) No. 3, 169-201; 19616, Astron. J. 66, 514-516; 1962, J. Res. NBS, 66B (Math, and 

Math. Phys.) No. 1, 5-13. 

(Paper 67B4-103) 



222 



